Вестник ДГТУ. 2012. № 7 (68) 








УДК 539.376 
Модели теории ползучести бетона и их конечноэлементная реализация 
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Разработан конечноэлементный метод решения плоской задачи теории ползучести бетона с учётом старе- 
ния. В качестве наследственных функций второго рода использованы выражения, базирующиеся на форму- 
лах для меры ползучести, предложенных Н. Х. Арутюняном и С. В. Александровским. С помощью символьно- 
го процессора получены выражения для ядер релаксации бетона, удобные для программирования. Разрабо- 
таны устойчивый шаговый алгоритм и соответствующая программа, позволяющие выполнять расчёты при 
переменном нагружении с учётом быстро набегающей ползучести в начальный момент нагружения и ча- 
стичной обратимости деформаций ползучести при разгрузке. На основе имеющихся экспериментальных дан- 
ных по ползучести центрально сжатых призматических бетонных стержней и изгибаемых железобетонных 
балок осуществлена верификация разработанного математического и программного обеспечения. 
Ключевые слова: метод конечных элементов, плоская задача, наследственная теория старения, численное 
интегрирование. 


Введение. В настоящее время накоплен значительный экспериментальный и теоретический ма- 
териал по ползучести бетона [1, 2, 3]. Вместе с тем известные подходы к расчёту бетонных и же- 
лезобетонных конструкций с учётом ползучести ориентированы главным образом на решение за- 
дач с относительно простой геометрией изделия [1, 3, 4]. В работе [5] приведён конечноэлемент- 
ный алгоритм решения задачи теории упругой наследственности. Однако данная теория постули- 
рует полную обратимость деформаций ползучести, что ограничивает её применение рамками 
только «старого» бетона. Поэтому актуальной является проблема разработки более общего ко- 
нечноэлементного алгоритма, реализующего модель теории наследственного старения, позволя- 
ющую моделировать процесс непрерывного бетонирования, а также рассчитывать потерю пред- 
варительного натяжения арматуры, обусловленную ползучестью стареющего бетона. 
Наследственные функции второго рода для бетона. Физические соотношения для плоской 
задачи наследственной теории старения запишем в матрично-операторной форме: 


{6} = [Е (#)1-В) {Е} 


где векторы-столбцы напряжений и деформаций {0} ={010,0,}, (= {Е =,2=,} (Т — 
[4 
символ транспонирования); [Е (Е) — матрица упругости материала; В=, = [в (Е, т)=, (т)ат, 


К (Е,т) — ядро релаксации (наследственная функция П рода). Здесь введены обозначения: 


Т — «возраст» материала в момент приложения нагрузки; Е — временная координата, от- 
считываемая от Т до некоторого текущего момента времени наблюдения. Параметры ЕЁ и т 
измеряются в сутках. 

Вид функции К (Ё,т) базируется на принятой механико-математической модели ползуче- 


сти бетона и применяемой функции меры ползучести № т В рамках теории старения выраже- 


ние для ядра релаксации К (Ё‚Т) имеет вид [1, 4]: 


99 


Технические науки 








где функция меры ползучести, предложенная Н. Х. Арутюняном [2], 


М(Е, = [< а е"* ый 


модуль упругости материала Е (#) =Е, пе": С, А,, \, Е, В — константы, определяемые 


из опытов на ползучесть при одноосной деформации. В формуле (1) величину т, можно считать 


равной моменту времени распалубки. Достоинством представления наследственной функции в 
форме (1) является то, что после обработки этого выражения в среде символьного процессора 
компьютерной математики Маре получаем следующую легко программируемую формулу: 


В (Е т) = 1 о и = 


_Е() 
Е (т) [ХзХовэ, (# ‚тв, (Е т)Ву =, (т)у(В+у) +, (Е, _. (2) 
+Х.Х. ©, (т) |, (Е) (В+у)- в, (Е, т)у-ц, (т) (В+у) +, о 
хХ,Х, в, (Ё,т) |, (Е) (В+ у) в, (Е, т)у —ц, (т) (В+) +6, ‚т)у], 


где введены обозначения: 





Е’(т) = Е Ве" ; х, =е ие Х› = Е, (т,С, +А,); те 
у. (х) = еУ(**",) ; (Ё,т) = е\*+т,)-вт ; (Е,т) с ет) 3 (Е, т)=е -У(т+6) , 


Графики функции К(Ет) для бетона, построенные с помощью м (2), для различных зна- 


чений параметра т приведены на рис. 1. На этом рисунке линиям 1, 2, 3, 4, 5, 6 соответствует 
«возраст» бетона (в сутках) 2, 8, 14, 20, 30, 40. Значения констант принимались равными [1]: 


№ =9,9388-10—11 (Н/м?) 1; д, = 47095-10! сут./(Н/м?); Е, = 2,55-1018 Н/м?; у = 0,03 сут. "; 


В = 0,206 сут." 

Всесторонний анализ недостатков теории старения, главным из которых является полное 
отрицание обратимости деформаций ползучести, представлен в работе [1]. Следующим шагом 
развития теории ползучести бетона является модель упруго-ползучего тела, или теория наслед- 
ственного старения. Суть теории наследственного старения состоит в предположении о частичной 
обратимости деформаций ползучести при условии выполнения принципа наложения воздействий. 
Наиболее точно экспериментальные данные о ползучести бетона описываются с помощью выра- 
жения для меры ползучести, предложенного С. В. Александровским [1]: 


С(Ет)=ч(т)- [5 А. А(т)1 - г" , (3) 


А А.- 
Здесь обозначено: ЕС А(т)=С,-С, + А . Величины А,, А,, А,, С,, С., \, 








а — опытные константы. 
Выражение для наследственной функции К(Е,т), полученное на основе формулы (3), 


имеет вид [1]: 
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К (Е,т) = =! (т)-'(т)(е\ -^,)- К’ (т)- [К (т) (е" -А,)е"®] й 
(4) 


Е 


х| К (т)='(т)е" Чат + В, (ем 











пе ется КО торе! "9-ГО" Ат 
(#2) (е" [ЕК а 
(= о Оч" А) [2 @)-©(@)- 
А (Е 





-аЕ (©) [Е (#)А(#)]| - @Е? (Е) (: ) не (Е )(е* - д, )[Е? (#)-К?(#)] + 
к(Е) 2 


+2аЕ? (Е) А (К) Е’(Е)(е" -А,}. 
После обработки выражения (4) символьным процессором системы Маре и группировки 
членов получим следующую формулу для наследственной функции К (Е,т) : 


кол) (т) ='(т)(е" сек ие. 


{Ее м Ё суе- АЕ м + (ТЕ (т) + (5) 





т 


Е (т)уе"е" Ио е-@\е-т) 
] (98.0 , 


"1+ А(Т)Е(т 






































й 4 2А, ь 2А.уе" (6)у’е” , 23, (ЕК)у’е”“ и. ы 
н(Е) В. (Е) [а в(е*-А) (е" А.) (е"- А, | А, 
{7-20 пе ее") ЕЕ -= ЕО, 
роны Е.В В н- 

Е (А д) 2(А, да Е) В?е пе 9 а 
ЕВЕ ЕО (Е) (е* ный -К? (Е 

а (#) АЕ (-(а-^) Е НЯ ем А,)х 

ЕЕ) зу (Е) 3 Аа 
х2| Е (Е)Е ИО ты ( Е (К л(Е)Е’ (6 (е -^) 
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А, - А, Е (т) А 


9, (т) =1+А(т)Е (т); 3, (т)=А(т)Е'(т) 5 ; 9. (т)=С,+—. 





= 

Интегралы п(т) и и(т), входящие в выражение (5), вычисляем численно с помощью 
формулы трапеций. 

Графики функции К (Ё,‚т) для бетона различного «возраста», построенные на базе фор- 
мулы (5), показаны на рис. 2. Здесь номера линий соответствуют значениям параметра т (в сут- 
ках): 1 — т=2;2 — т=4; 3 — т=6б; 4 — т=10; 5 — т=20; 6 — т=30. Значения констант 
принимались равными ПЕ № = 9,9388-107"" (Н/м?)"; №, = 7,7064-107" (Н/м?)"; 
А, = 4,7095-10`1° сут./(Н/м”); А, =1; А, =3,4822.10— сут./(Н/м?); Е‚ =2,55:100 Н/м?; а=6 сут. "; 
У = 0,03 сут. "; В = 0,206 сут. 1 


И 








[ т т т т | т ] Ц т 
0 20 `60 100 '140 `180; сут 00 20 3% 40 
Рис. 1. Графики функции К (Е т) в теории старения Рис. 2. Графики функции к(Ет) в теории наследственного 
старения 


Отличительная особенность графиков рис. 2 от графиков рис. 1 — резкое убывание функ- 
ЦИИ К (Е,т) в момент времени Ё=т, что с физической точки зрения отражает наблюдаемое на 
практике явление «мгновенной» ползучести, условно причисляемое к упруго-мгновенной дефор- 
мации в момент загружения. 

Для описания ползучести «старого» бетона можно использовать теорию упругой наслед- 
ственности с ядром релаксации типа [1]: 


Е (Е, т) = Ае® 7 + Ае 9, 


Е 
где А, Е + ал, )[Е› (уф+ ад, ) фр» | +а?А, + уф]; А, = Еь (Уфо +аАь) — А; 


1 2 





1 
р: .= 5 +у+ Е, (УШь + ад, ) + ЧЕ? (Уфо +ал, + (а -5) —2Е‹ (Уф, +ал, )(а-у).; 


Фо =Л; А, =С, - С, . 
Конечноэлементная реализация плоской задачи теории ползучести. Линейно-упругие 
перемещения в глобальной Декартовой системе осей {27,,7,} в произвольной точке конечного 


р о т Е 
элемента (КЭ) задаём с помощью вектора-столбца перемещений {и} = {ии} . Введём векторы- 
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столбцы наследственных деформаций {&}, перемещений {0} и узловых перемещений {и’}: 


{Е} =(1-в){=}, {2} =(1-в) {и}, {и} =(1-в) и}. 
Векторы-столбцы {&} и {@} имеют структуру, аналогичную {=} и {и}. Размерность векто- 


ра-столбца {и/} определяется числом узлов КЭ. Установим связь между {й} и {и/} в виде 


{1} =(1- в [ЕР {и}, 
где [2] — матрица, образованная из функций формы КЭ. Отметим, что для плоского восьмиузло- 


вого КЭ матрица [Е | имеет размерность 2х8. 


В соответствии с принципом возможных перемещений получим операторно-матричное 
уравнение 
[«а-в){и}- {и} =0, (6) 


где матрица жёсткости | К | и вектор-столбец узловых сил {/} КЭ 
[е= [2 (ЕС, (= ГДЕ ада + [27 {р}; 


т уз ВН К 5 К 
{9} ={9,9›} , {р} ={р.р›} — векторы-столбцы объёмной и распределённой нагрузки, задавае- 
мой в глобальном базисе; и, — объём, занимаемый КЭ; $. — поверхность КЭ, к которой прило- 
жена распределённая нагрузка. Матрица [2] для восьмиузлового КЭ имеет размерность 3х16. 


Выражение [2] для плоского полилинейного КЭ получено в работе [6]. 


[4 
Для вычисления интеграла В {и} = [К(Е,т) {и (т)} т воспользуемся численным методом, 


1 


основанным на формуле трапеций. Разобьём рассматриваемый временной интервал 77] на т 


равноотстоящих временных шагов ДЁ так, чтобы Ё = ТДЁ. Тогда выражение (6) можно записать 
в форме 
Ви =в (Е.Е) м} ДЕ/2+ Ув(Е (т - ЛА Е+в(Ет) {и} ДЕ /2-{"}=0 
7= 


или в компактном виде 
[К и} = {и} + [® Ци} + [К, ] и, ве ‚ (7 

где [к„|=[]а-в(.8)4Е/2; [К] =[К] а - ВЕ, т) /2; [к,|=[ К] а-в(Е(т- ЛАВ Ае. 
В выражении (7) вектор-столбец {и} соответствует упруго-мгновенному решению зада- 


чи. Рассмотренная шаговая конечноэлементная процедура реализована на языке ГОКТКАМ. В 
расчётных моделях полилинейные КЭ можно ансамблировать с плоскими стержневыми КЭ балоч- 
ного типа. 

Числовые примеры. Мгновенное нагружение и последующая ступенчатая разгрузка призмати- 
ческого бетонного образца с размерами 6х6х30 см. Характеристики материала были приведены 
ранее. «Возраст» бетона в момент нагружения образца — т = 28 сут., момент времени распалуб- 
ки — т, = 2 сут. Учитывая симметрию задачи, рассмотрим > часть образца и применим равно- 


мерную конечноэлементную разбивку 2х8 (2 КЭ по ширине и 8 КЭ по высоте). Шаг интегрирова- 
ния ДЁ принимаем равным 1 сут. График полного перемещения нагруженного торца образца на 
интервале наблюдения Е = 28...200 сут. приведён на рис. 3. 
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20 60 100 140 180 г, сут. 


Рис. 3. Ступенчато-убывающий закон нагружения призматического образца: 
1 — теория упругой наследственности; 2 — теория старения; 3 — теория наследственного старения 


Как видно из представленных данных, для модели упругой наследственности (линия 1) 
выполняется постулат о полной обратимости деформаций ползучести. Для модели упруго- 
ползучего тела (линия 3) характерна частичная обратимость неупругих деформаций ползучести, 
что соответствует результатам аналогичного натурного эксперимента А. Д. Росса [1, 3]. 

Ползучесть однопролётной железобетонной балки, мгновенно нагруженной сосредоточен- 
ной силой посередине пролёта (рис. 4). Размеры балки приведены в метрах. Диаметр арматуры — 
8 мм. Модуль упругости арматуры 2,1-10° МПа. Разбивку балки на плоские КЭ выполняем сеткой 
6х40 КЭ (6 КЭ по высоте и 40 КЭ по длине). Арматуру моделируем стержневыми КЭ (40 КЭ). Шаг 
интегрирования — ДЁ = 1 сут. 


Р=3,6 кН 





Рис. 4. Расчётная схема однопролётной железобетонной балки 


Рассмотрим два варианта армирования: Г — арматура расположена в нижней (растянутой) 
половине балки, П — симметричное армирование растянутой и сжатой частей балки. Результаты 
численного моделирования процесса ползучести в виде графиков зависимости прогиба Ё› в цен- 


тре пролёта от времени наблюдения Ё и «возраста» бетона т представлены на рис. 5. 

На рис. 5 штриховые линии относятся к варианту армирования Т, сплошные линии — к ва- 
рианту армирования 11; номера линий 1 и 2 соответствуют значениям параметра т, равным 14 и 
28 сут. Момент времени распалубки т, = 2 сут. Значения упруго-мгновенного прогиба (т) для 


вариантов армирования 1 и П соответственно составляют 0,4890-10`3 м и 0,4199-103 м. 

Как видно из представленных на рис. 5 графиков, наличие арматуры в сжатой зоне балки 
обусловливает снижение ползучести для всех реализованных моделей бетона. Вместе с тем толь- 
ко модели теории старения (рис. 5, 6) и теории наследственного старения (рис. 5, в) описывают 
наблюдаемую на практике быстро натекающую деформацию ползучести /,, (т+1), возникающую 


тах 
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сразу после упруго-мгновенной деформации балки. Значения #,, (#)/*.., (т) (#)/Е„„ (т+1), 


где { = 150 сут. для теории старения и теории наследственного старения приведены в =. 


т 








0 50 100 150 0 50 100 150 
т, сут. 1, сут. 


0 
—0,6 


. 
ВА 
ЕР 
>< 
\ 


Е ° 10 


0 50 100 150 
р, сут. 
в) 
Рис. 5. Графики /„„(Ё) для железобетонной балки: 


а — теория упругой наследственности; 6 — теория старения; в — теория наследственного старения 

















Значения Ё,, (Е (Е)/Ё.. (ТИ, (Е)/1.„ „ (Т+1) для однопролётной балки 
Е. (ЕЕ, (Т (т + 1) 
рр и.) =. (Е)/"... 
т = 14 су. Т = 28 сут. т = 14 су. Т = 28 сут. 
Теория старения 3,19 3,09 1,494 1,494 
Теория наследственного 3,87 3.21 1,63 1,39 
старения 

















Из приведённых в табл. данных следует, что значения РЁ 


ий и  (Т+1), соответствую- 
щие теории старения, для моментов времени т = 14 сут. ит = 28 сут. совпадают. Тогда как в 
действительности по мере «старения» бетона это отношение должно уменьшаться. Отметим, что 
данный эффект адекватно описывает теория наследственного старения. 

Выводы. 

1. В рамках теории старения и теории наследственного старения с помощью символьного 
процессора компьютерной математики Маре получены выражения для ядер релаксации, удобные 
для программирования. 

2. Разработана и программно реализована шаговая процедура метода конечных элемен- 
тов, позволяющая моделировать процессы последействия в бетонных и железобетонных кон- 
струкциях с учётом старения бетона. 
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3. На тестовых примерах выполнена численная апробация разработанной конечноэле- 
ментной программы. 
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СОМСКЕТЕ СВЕЕР ТНЕОКУ МОБЕТ$ АМО ТНЕТВ ЕТШТЕ- 
ЕГЕМЕМТ ТМРЕЕМЕМТАТТОМ 


Р. Р. Сау4тВигом 

(Роп За Тесптка! Упмег$КУ), 

Е. В. 15КВаКома 

(боиё-Кизап $(ае Тесптса! Упмег$Ку) 


Тре Ипйе-е/етепЕ тето4 Гог йе р/апе рго Мет 5о/ийоп о! {пе сопсгее сгеер {теогу ий гедага ю ада 1/5 аЕ- 
уебреа. Тйе ехргез$!оп5 Базеа оп те югттиа$ Гог сгеер уа/иез ргорозеа Бу М. К. Апфуипуап апа $. И АЕхап- 
СГоу$Ку аге изе4 аз те Вегеаагу ГипсНоп$ ог {Те 5$есопа Кта. Тпе Репа ргодгатт/пд ехргеззопз юг {те соп- 
сгее ге/ахавоп Кегпеб аге оматеа Бу те утро/с ргосе5зог. Тве аигаМе %ер-Бу-$ер а/догЁйт апа 1е аррго- 
рпаёе ойи/аге Юг реготттпа ассоипйпд ипдег пе ата Ме !оа@тд ий" гедага №0 {те ГазЁ тодепЕ сгеер аЁ пе 
та! юаа апа рагИа! геиегЮИЙу ог {пе сгеер Пои/ ипаег аЕоаатд аге аеуеоред. Тпе уетсайоп оЁ {пе дЕуе!- 
ореа зоймаге 15 сагтеа оиЁ оп {пе Базе о! те аиаЙа ее ехрептега! ааа оп #1е сгеер оЁ те !есйу сотргезме 
ризтайс сопсгее гоа5 апа Бепатд сопсгее Беатс. 

Кеуигога&: ЕпПе-е/етепЕ тетод, рапе рго Мет, Вегеайагу адта #пеогу, питепйса! п{едгайоп. 
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